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THE  ACCURACY  OF  THE  SUBSTITUTE-STRINGER  APPROACH  FOR 
DETERMINING  THE  BENDING  FREQUENCIES  OF 
MULTISTRINGER  BOX  BEAMS 
By  William  W.  Davenport 

SUMMARY 


The  accuracy  of  the  substitute-stringer  approach  for  including  the 
effects  of  shear  lag  in  the  calculation  of  the  transverse  modes  and  fre¬ 
quencies  of  multistringer  box  beams  is  investigated.  Box  beams,  the  covers 
of  which  consist  of  nomal-stress -carrying  stringers  on  sheets  which  carry 
not  only  shear  but  also  normal  stress,  are  analyzed  exactly.  Frequencies 
of  beams  with  various  numbers  of  stringers,  obtained  by  means  of  this 
exact  analysis,  serve  to  determine  the  possible  accuracy  of  the  frequencies 
obtained  by  the  substitute-stringer  approach . 


INTRODUCTION 


The  use  of  a  substitute-stringer  approach  for  including  the  effects 
of  shear  lag  in  the  vibrational  analysis  of  built-up  box  beams  was  inves¬ 
tigated  in  reference  1.  Various  thin-walled  rectangular  tubes  were  ideal¬ 
ized  to  substitute-stringer  structures  and  the  frequencies  of  the  first 
few  bending  modes  of  the  idealized  structures  were  compared  with  those 
of  the  original  beams.  The  results  indicated  where  the  substitute  stringer 
should  be  located  in  order  to  take  into  account  accurately  the  effect  of 
shear  lag  on  the  frequency;  for  a  fairly  wide  variety  of  tube  proportions, 
locating  the  substitute  stringer  midway  between  the  web  and  the  centroid 
of  the  half -cover  yielded  frequencies  accurate  within  1  or  2  percent  for 
higher  as  well  as  lower  modes  of  vibration.  However,  because  the  rec¬ 
tangular  tubes  employed  in  the  investigation  had  constant  wall  thickness 
there  was  no  definite  information  given  regarding  the  accuracy  of  fre¬ 
quencies  that  may  be  obtained  for  a  beam  which  has  discrete  flanges  and 
stringers  in  the  covers  and  different  web  and  cover-sheet  thickness. 

The  purpose  of  this  report  is  to  examine  the  accuracy  of  the 
substitute-stringer  approach  when  applied  to  a  beam  more  nearly  repre¬ 
sentative  of  actual  built-up  box  beams.  Frequencies  of  various  realistic 
multistringer  beams  of  the  type  shown  in  figure  1(a)  are  determined  by 
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means  of  the  substitute-stringer  approach  as  given  in  reference  1.  These 
frequencies  are  then  compared  with  the  frequencies  obtained  by  means  of 
an  exact  vibrational  analysis  of  a  multistringer  beam.  Conclusions  are 
thereby  made  with  regal'd  to  the  proper  location  of  the  substitute  stringer. 


SYMBOLS 


A 

Ap 

al 

AO 

AP 

AS 


a 

b 

b0  =  a 
bp 

^r+l 


*>S 

*>C 

Bi 

C 

Ci 

Di 


(4) 


(4) 

(*) 


total  cross-sectional  area  of  stringers  in  a  half -cover  of 
multi stringer  box  beam 

cross-sectional  area  of  flange  of  substitute-stringer  structure 
cross-sectional  area  of  a  substitute  stringer 
cross-sectional  area  of  flange  of  multistringer  box  beam 
cross-sectional  area  of  pth  stringer  of  multistringer  box  beam 
effective  shear-carrying  area 

half -depth  of  beam 
half -width  of  beam 

width  of  pth  panel 

half -width  of  middle  panel  when  N  is  even,  nonzero  constant 
when  N  is  odd 

distance  between  web  and  adjacent  substitute  stringer 
distance  between  web  and  centroid  of  area  of  half -cover 

parameter  defined  by  equation  (A17)  or  equation  (MO) 
constant 

parameter  defined  by  equation  (Al8)  or  equation  (Ml) 
parameter  defined  by  equation  (A25)  or  equation  (M2) 


3656 


modulus  of  elasticity 


1 

OtcX,2^  +  (r+l),q) 

shear  modulus  of  elasticity  (taken  equal  to  E/2.65  herein) 


bending  moment  of  inertia 


♦ 5 


Oj)  +  ^ 


frequency  coefficient,  o>U==- 


1  Ittt 

coefficient  of  shear  rigidity,  -  W — 


s)2  .  kBV 


half-length  of  beam 

number  of  stringers  in  a  cover 

parameter  defined  by  equations  (A3^)  ®®d  (A55)  or  t>y  equa¬ 
tions  (AM-3)  and  (AMO 


determinant  given  by  equation  (A30) 


t(j  cover-sheet  thickness 

tc'  effective  cover-sheet  thickness  for  normal  stress 

ty  web  thickness 

T  maximum  kinetic  energy 

U  maximum  strain  energy 

Up  longitudinal  displacement  of  a  point  on  the  pth  panel 

w  vertical  displacement  of  cross  section  of  beam 

x  longitudinal  coordinate 

yp  chordwise  coordinate  for  pth  panel 

Fourier  series  coefficients 


IjJjPjq.;111*11  integers 


B_  ■  — 3~ 

*  2GtcL 

6.^  Kronecker  delta  (0  when  i  ^  j;  1  when  i  »  j) 


M{p) 

n 

Ui 

cue 


Lagrangian  multiplier 
mass  of  beam  per'  unit  length 

parameter  defined  by  equation  (A8) 
natural  frequency 

natural  frequency  of  multistringer  beam  in  figure  3 
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THE  MULTISTRINGH*  BOX  BEAM  AND  THE  SUB8TITUTB-8TRINGER  STRUCTURES 


A  doubly  symmetrical  multistringer  box  beam  of  the  type  considered 
in  the  present  paper  is  shown  in  figure  l(a).  In  order  that  the  covers 
of  this  structure  may  behave  realistically  with  regard  to  shear-lag 
effects,  the  cover  sheets  are  permitted  to  carry  not  only  shear  but  also 
normal  stress.  The  following  simplifying  assumptions  are  made: 

(1)  The  flanges  and  stringers  carry  only  noxrnal  stress. 

(2)  The  webs  carry  only  shear  (the  bending  resistance  of  the  web 
is  included  in  the  flanges). 

(3)  Longitudinal  inertia  (rotary  inertia)  is  neglected. 

(4)  The  cross  sections  maintain  their  shape. 

Assumptions  (2)  and  (3)  are  generally  known  to  be  good  for  reasonably 
shallow  beams  such  as  those  considered  in  this  report;  assumption  (4) 
is  good  if  a  normal  amount  of  bulkhead  stiffness  is  present. 

In  the  appendix  an  exact  vibration  analysis,  similar  to  those  of 
references  1  and  2,  is  carried  out  for  a  multistringer  beam  of  the  type 
shown  in  figure  l(a).  The  Rayleigh-Ritz  energy  procedure  is  used  in 
conjunction  with  appropriate  Fourier  series  and  Lagrangian  multipliers 
to  obtain  frequency  equations  for  a  box  beam  with  any  number  of  stringers. 
Although  the  analysis  in  the  appendix  allows  unequal  spacing  of  the 
stringers,  only  the  case  of  equally  spaced  stringers  is  considered 
hereinafter . 

The  substitute-stringer  structure  for  the  beam  of  figure  1(a)  is 
shown  in  figure  l(b).  lhe  flange-web  combination  of  this  structure  is 
the  same  as  that  of  its  prototype;  the  covers,  however,  consist  of  sub¬ 
stitute  stringers  which  carry  only  normal  stress  and  sheets  which  carry 
only  shear.  The  magnitude  of  the  shear-lag  effect  in  this  structure 
depends  on  the  location  of  the  substitute  stringers.  A  suggested  stringer 
location  given  in  reference  1  is 

bg  -  0.5bc  (l) 

where  bg  is  the  distance  between  a  flange  and  the  adjacent  substitute 
stringer  and  bc  is  the  distance  between  a  web  and  the  centroid  of 
normal-stress-carrying  area  of  the  half-cover. 

A  vibration  analysis  of  the  substitute-stringer  structure  is  carried 
out  in  appendix  B  of  reference  1. 
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COMPARISON  CP  PRBCfJHJCIIS  AS  OBTAINED  BY  THE  RJIiTIBTKWaR 
SOLUTION  AND  THE  SUBSTTIUTE-STRINQER  APPROACH 


In  order  to  investigate  the  possible  accuracy  of  the  substitute- 
stringer  approach  for  determining  the  frequencies  of  multistringer  box 
beams,  comparisons  are  made  between  the  frequencies  of  four  multistrlnger 
box  beams  as  obtained  by  the  analysis  given  in  the  appendix  and  the  corre¬ 
sponding  frequencies  of  their  substitute-stringer  structures  as  obtained 
by  the  analysis  given  in  reference  1.  The  first  three  symmetrical  trans¬ 
verse  modes  of  each  beam  with  free -free  end  conditions  are  considered. 

The  four  beams  have  cross  sections  as  shown  in  figure  2  and  are 
identical  except  for  the  manner  in  which  the  total  stringer  area  is 
distributed.  Proportions  caramon  to  all  four  beams  are 


tyr 

-T  *  1-25 

tC 


where  L,  b,  a,  ty,  and  t^  are,  respectively,  the  half-length,  the 

half -width,  the  half -depth,  the  veb  thickness,  and  the  cover-sheet  thick¬ 
ness;  the  total  cross-sectional  area  of  the  stringers  in  a  half -cover  is 
given  by  A  and  the  cross-sectional  area  of  a  flange  is  given  by  Aq. 

In  addition  to  the  preceding  properties,  the  four  beams  are  alike  in  that 
stringers  of  each  beam  are  equally  spaced  across  the  covers.  The  beams 
thus  differ  only  in  the  number  of  stringers  on  a  cover;  the  cases  con¬ 
sidered  are  for  N  =  1,  4,  7>  and  «. 

In  figure  2  the  various  relations  needed  to  obtain  the  parameters 
of  the  multistringer  solution  for  the  four  beams  are  given.  When  N  is 
infinite  the  stringers  are  smeared  out  to  a  "stringer  sheet"  for  which 
the  effective  normal-stress -carrying  thickness  t^'  is  different  from 

the  actual  cover-sheet  thickness  tg. 

The  cross  section  of  the  substitute-stringer  idealization  is  also 
shown  in  figure  2.  In  accordance  with  the  idealizing  procedure  of  the 
substitute-stringer  approach,  the  following  relations  are  used  to  deter¬ 
mine  the  parameters  needed  in  the  substitute-stringer  solution  of  refer¬ 
ence  1: 

Ay  *  Aq  Ajj  -  A  +  bt(j 

where  Ay  and  A*,  are,  respectively,  the  cross-sectional  area  of  a 

flange  and  stringer  of  the  substitute-stringer  beam.  It  should  be  noted 
that  the  only  geometric  property  of  the  substitute-stringer  structure 
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which  may  he  different  for  each  of  the  four  beams  is  the  value  of  bg. 

For  the  substitute -stringer  location  defined  by  equation  (l),  bg  varies 
by  virtue  of  the  variation  of  bg.  The  values  of  be  for  the  four  beams 
are: 


For  1  stringer, 


For  4  stringers, 


For  7  stringers, 


For  infinity  stringers, 


bc  *  0.714b 
be  *  0.545b 
bc  *  0.551b 


bc  =  0.500b 


The  frequency  coefficients  kj$ 
substitute-stringer  solution  where 


obtained  from  both  the 


bg  *  0.5bc  and  the  multi  stringer 


solution  for  each  of  the  four  beams  are  given  in  table  I.  In  addition 
to  these  results,  the  frequency  coefficients  obtained  by  a  solution 
which  includes  transverse  shear  but  not  shear  lag  have  been  included. 
This  latter  set  has  been  presented  in  order  to  demonstrate  the  influence 
of  shear  lag  and  was  obtained  from  equation  (15)  of  reference  5  with  the 
rotary  inertia  parameter  kgj  equal  to  zero. 


The  influence  of  stringer  location  on  the  frequency  co  of  the 
substitute-stringer  structure  for  each  of  the  four  beams  is  presented 
graphically  in  figure  5-  The  frequency  is  given  in  the  form  of  its 

relative  error  -  1  when  compared  with  the  exact  frequency  of 
(be 

the  multistringer  structure;  the  curves  indicate  how  the  errors  vary 
with  substitute-stringer  location. 


DISCUSSION  OF  RESULTS 

The  results  given  in  table  I  show  good  agreement  between  the  fre¬ 
quency  coefficients  as  obtained  by  the  substitute-stringer  approach  with 
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bg/bc  *  0.5  and  those  obtained  by  the  multistringer  solution.  However, 
as  may  be  noted,  the  frequency  coefficients  for  the  substitute  structure 
with  bg/bc  *0.5  are  slightly  high  for  all  but  the  third  modes  of  the 

cases  for  N  ■  1  and  N  *  ».  Reference  to  figure  5  indicates  that,  for 
all  but  the  case  of  N  *  *,  the  maximum  errors  for  the  modes  considered 
would  be  reduced  by  using  a  value  of  bg/bc  *  0.55>  for  example,  for  the 
case  of  N  =  7,  the  maximum  error  for  this  new  value  would  be  less  than 
1  percent. 

The  results  of  table  I  indicate  that  the  stringer  area  distribution 
has  little  influence  on  the  frequency;  however,  the  assumption  should 
not  be  made  that  the  influence  of  shear  lag  is  negligible.  As  may  be 
seen  in  table  I,  the  reduction  in  frequency  due  to  the  inclusion  of  shear- 
lag  effects  is  4.8,  16.7,  and  19.1  percent,  respectively,  for  the  first 
three  symmetrical  modes  of  the  case  where  N  =  7* 

CONCLUDING  REMARKS 


The  numerical  results  of  the  present  paper  indicate  that  the 
substitute-stringer  approach  can  yield  accurate  frequencies  for  multi¬ 
stringer  box  beams.  Review  of  the  numerical  results  of  NACA  Technical 
Note  3158  together  with  those  of  the  present  paper  suggests  that  a 
value  of  O.55  for  bg/bc  (ratio  of  the  distance  between  web  and  adjacent 

substitute  stringer  to  distance  between  web  and  centroid  of  the  area  of 
the  half-cover)  defines  a  slightly  more  appropriate  substitute-stringer 
location  for  built-up  box  beams  than  bg/bc  =  0.50,  the  value  suggested 

in  NACA  Technical  Note  3158- 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va. ,  January  6,  1956. 
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APPENDIX 

VIBRATION  SOLUTION  OF  A  MULTISTRZNGER  BOX  BEAM 


The  simplifying  assumptions  for  the  multistringer  beam  are  given 
in  the  body  of  the  paper  and  a  cross  section  consisting  of  two  cases  is 
shown  in  figure  4.  Case  I  applies  vhen  N  is  even  and.  case  II  vhen 
N  is  odd,  where  N  is  the  number  of  stringers  in  a  cover.  In  addition 
to  the  dimensions  and  coordinates  defined  by  figure  4,  the  length  of  the 
beam  is  assigned  the  value  2L,  and  x  is  defined  as  the  longitudinal 
coordinate  with  its  origin  at  the  midpoint  of  the  beam.  It  should  be 
noted  that  br+i,  the  distance  between  the  rth  stringer  and  the  center 
line  for  case  I,  is  half  the  width  of  the  middle  panel. 

For  a  transverse  mode  of  vibration  let  w(x)  be  the  amplitude  of 
a  vertical  displacement  of  a  cross  section  and  let  up(x,ypj  be  the 

amplitude  of  a  longitudinal  displacement  of  a  point  on  the  pth  panel. 
Then  the  maximum  strain  energy  is 


r+1 


U  =  2E 


Vl 


-f 

',P  J-I 


p=l  1  +  ^(r+lbp^-iA  dx 


-  8s6(r+l P  ^  ^  + 


2Gtr 


-r  (s  ■  ^ « 


(Al) 


where 


»,  ■ 1  ■- 
W  2 


(A 2) 


and 


6(r+l),p  *  0 
6(r+l),(r+l)  “  1 


(p  /  r  +  1)  (A3) 
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The  maximum  kinetic  energy  when  the  influence  of  longitudinal  inertia 
is  neglected  is 


T  =  i  /  w%x  (A4) 

2  J-L 

where  p  is  the  mass  per  unit  length  of  the  structure  which  the  multi¬ 
stringer  beam  represents.  It  should  be  noted  that  equations  (Al)  to  (Ah- ) 
hold  for  N  =  0. 

In  the  following  Sections  appropriate  infinite  trigonometric  series 
are  assumed  to  represent  the  displacements  and  then  the  Rayleigh-Ritz 
energy  procedure  is  applied. 


Symmetrical  Transverse  Modes  of  a  Free-Free  Multistringer  Beam 

Appropriate  trigonometric  series  for  the  symmetrical  transverse 
modes  of  a  free-free  beam  are 


up(x,yp)  -  Z  %n(p)si”  g*  “8  T? 

m=l,3,5  n=0,l,2 


2L  -  bp 

(p  =  1,2,  .  .  .  (r  +  1)) 


(A5) 


w(x)  =  C  +  z.  cm  cos  Sp  (A 6) 

m-1,3,5 


The  choice  of  these  particular  series  was  guided  by  the  orthogonality 
required  for  simplification  of  the  expressions  in  the  strain-energy 
equation;  C  in  equation  (A 6)  is  introduced  in  order  to  allow  w(±L) 
to  be  unrestricted. 


Because  of  the  geometry  of  the  structure  the  following  restraining 
relations  must  hold: 

Up(x,bp)  -  up+L(x,0)  =0  (p  ■  1,2,  .  .  .  r)  (A7) 

Substituting  equation  (A5)  into  equation  (A7),  multiplying  by  sin  gj*- 
where  i  *  1,3,5>  •  •  •  >  and  integrating  fran  -L  to  L  yields 

Z  Ei„(p)(-Un  -  Hn(P+13  -  *l(P)  - 0 

n*0,l,2L" 

(p  “  .  .  .  r)(i  ■  1>3>5>  •  •  •)  (A8) 
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In  accordance  with  the  Rayleigh-Ritz  procedure  the  expression 

u  -  T  -  ±  t.  >»(P)fe(P)  (A9) 

p=l,2  m*l,3,5 


is  held  stationary  with  respect  to  the  coefficients  of  the  assumed  series. 
The  A's  in  expression  (A9)  are  Lagrangian  multipliers  introduced  to 
maintain  the  condition  stated  in  (a8).  The  expressions  for  U  and  T 
of  equation  (A9)  are  given  by  equations  (Al)  and  (A4),  respectively,  and 
the  displacements  are  in  turn  given  by  equations  (A5)  and  (A6).  Thus, 

differentiating  expression  (A9)  with  respect  to  the  ai^'s,  ci's, 
and  C  independently  and  setting  the  respective  results  equal  to  zero 
yields  the  following  equations: 


Klj(l)aij(l)  +  2Ri(l) 


Si  Mn(l)  +  fr  H(^)ci  -  =  0 

n-0,1,2  L  '2/ 


(i  =  1,3,5,  •  •  *)(J  —  0,^,2,  .  .  .) 


(A10) 


Si  HnM  ‘ 

n=0,l,2 

Sq^l^t-l)^!  -  8(r+i),q)  -  =  0 

(1  *  1,3,5,  •  •  •  )()  -  0,1,2,  .  .  . )(q  -  2,3,  .  .  .  r  +  l)(r  i  1)  (All) 


kiCi  + 


i¥) 


X  ain 

nO,l,2 


(1)  - 


2kB2kS2 


(i  -  1,3,5,  •  •  •)  (A12) 


where 


KlJCq)  “  Pq2(^)2(l  +  #oj)  +  (J*)2 


(A13) 


(Alka) 


x 
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_  2  Ebg 
•«?- 

(Alto) 

%(<l)  “  +  ®Xq 

(Alto) 

r  2  “q-l^q 

1  0tcL2(l  ♦  ^(r+xj^n) 

(Al4d) 

tca 

(Alto) 

2rk 

kB2  -  w  1 

B  El 

(Al4f) 

kg  =  ■El-g 

GAgl/ 

(Al4g) 

Ag  =  toty 

(Alto) 

e„  --5a- 

*  2GtCL 

(Aito) 

ki  *  [(I1)  “  toV] 

(Alto) 

Dividing  equation  (AlO)  "by  Kj_j  ^  and  summing  over 

J  *  0,1,2,  .  • 

gives,  on  solving  for 


£ 

J-0,1,2 


J-0,1,2 


(i  *  1,3,5,  •  •  •)  (A15) 
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Similarly,  for  equation  (All): 


jLi  a  -  8(r+i)q)-  ci^3q*: 

0s®  9 


(<1)0  ^(q-i) 


(i  =  1,3,5,  •  •  *)(q  ■  2,3,  •  .  .  r  +  1)  (Al6) 


where 


£ 


(-D‘ 


Bi 


M  =  _ J=0A,2  Kjj 


+  2Ri^q^  ^  - h 

KU 


j =0,1,2  - 


(q=l,2,  .  •  •  r)  (AI7) 


and 


Ci 


(q) 


£ 


J=0,l,2  Kjj 


(L  -  0#(r+l),q)  +  ^i 


(q) 


75T 


J=^7E,2  K±j 
(q  =  1,2,  .  .  .  r  +  1)  (Al8) 


In  equation  (AI7)  the  definition  of  Bj/r+^  has  been  neglected  because 

when  q  =  r  +  1  in  equation  (Al6)  the  entire  term  containing  Bi^r+^ 
drops  out. 

By  substituting  equations  (A15)  and  (Al6)  into  equations  (A10) 
and  (All)  the  following  equations  may  be  obtained: 


H,(l)  -  -V 

KljUJ 


[fRi(l)ci(1)  -  i]  +  Pi^i(1)  [(-i)J  -  2R1(1)b1(1^ 


(i  =  1,3,5,  •  •  *)(j  “  0,1,2,  •  .  .)  (A19) 


HACA  TN  5636 


'  1 

aij(q)  "  ~ W  +^(4)li-i)J  -  2Rt  (q)Bl(i)]l 

Kid 

(i  a  1/3,5,  •  •  *)(d  “  0,1,2,  •  .  .  )(q  *  2,3,  »  .  .  r)  (A20) 

Equations  (A19)  and  (A20),  when  multiplied  by  (-1)^  and  summed  over 
J  =  0  y\)<L)  •  •  •  j  y  i  0  Id 

2  *  -Bi(l)(f1)s(f)cl  +  Pl[fi(1)  +  Dl(1]])'i(1) 

(i  -  1,3,5.  •  •  •  )  (A21) 

and 

J=o7l,2  U  -J 

(i  =  1,3,5,  •  •  •  )(q  =  2,3,  •  •  •  r)  (A22) 

In  order  to  obtain  equations  in  the  form  of  (A2l)  and  (A22)  it  should  be 
noted  that 

I?1"’'1''1  -  3  Ju  ■  -»■“ 

(q  =  1,2,3,  •  •  •  r)  (A23) 


C  _1  PR  (4)b(0  f 
=^t,2  ZT7^7  1  j-£t. 


(-1)J  _  „  (<l)  .  «  (l) 

2^J-Cl  +Dl 

(q  =  1,2,3,  .  .  .  r)  (A24) 


where 


(q.  *  1,2,3,  .  .  .  r)  (A25) 

Substituting  equation  (A15)  into  equation  (A12)  results  in 

[ki  -  +  ^(f)Bi(1)xi(1)  -  a*V(£)<-i>“r  £ 

(i  =  1,3,5,  •  •  •)  (A«6) 

and  substituting  equations  (A2l)  and  (Al6)  into  equation  (A8),  when 
P  =  1,  gives 


(i  »  1,3,5,  •  •  •)  (A27) 

Substituting  equations  (Al6)  and  (A22)  into  equation  (A8)  yields 

+  Ci(p)  +  Dl(p) + ^  0l(p+ijv(p) . 

-£±i  Bi(p+1)Xi(p+1)(i  .  8(r+1)>(p+1))  -  0 

(i  *  1,3,5,  •  •  • )(p  “  2,3,  •  •  •  r)  (A28) 

By  Cramer's  rule,  equations  (A26)  to  (A28)  may  be  solved  for 
in  terms  of  C/P]_;  if  in  the  resulting  expression  the  determinants  in 
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By  expanding  the  determinant  of  (A30)  by  the  first  column  the  following 
recurrence  relation  may  be  obtained: 


where 


Si 


(r) 


(p  “  ^>2>3>  •  •  •  r  —  l)  (A31) 


ClM  ♦  Dl(r)  +  Cl<rtl> 

Dr 


and 


Si^r+1^  =  1 

Now,  by  dividing  the  numerator  and  denominator  of  the  right-hand  member 
of  equation  (A29)  by  and  then  using  equation  (A3l),  the  following 

equation  results: 


|ki  -  2H(f)2Ci(l^|  +  2E(f)2 

(b,<1>)V2> 

f  (1)  (1)  t2  (2)1  (2)  W  (2)\2  (3) 

ci  +  Di  +  ~  c±  s*  -  rW  )  Sl 
_  D1  J  °1 

(i  -  1,3,5,  :  (A32) 

Repetition  of  this  procedure  ultimately  yields  the  continued  fraction 


a 


where 
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permits  equation  (A33)  to  be  written  as 


2kn2k«2te)(-l)  2  C 


ki  - 


(*)2£ 


(i=lj3,5>  •  •  •)  (A36) 


Substituting  equation  (A36)  into  equation  (A13)  yields  the  frequency 
equation 


2  2  2 
kBJ  1  +  2kB  kg 


\  (m.)2  ist 

L  (f)^-Hr 

m-1,3,5  V 


O  (A37) 


The  following  closed  forms  for  the  summations  involved  in  B^'  , 

and  may  be  obtained  by  methods  similar  to  those  given  in 

reference  h- : 


J  =071,2  '  2(f)Pp  tanh  (*)Pp 


(p  =  1,2,  .  .  .  r  +  1)  (A38) 


t  = _ 1 _ 

J=°'1’2  I4nj  P  8(?)Pp  Slnh  (f)Tp 


(p  =  1,2,  .  .  .  r  +  1)  (A39) 


By  use  of  equations  (A38)  and  (A39)>  equations  (A17),  (Al8),  and  (A25) 
may  now  be  written  as 


P^  tanh 


Pt>  tanh 
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(p)  „  V 


tanh 


WPp  (f)?p  tmh  («)PP  +  En< 


(p  «  1,2,3,  .  .  .  r)  (M2) 

From  equations  (A34),  (A35),  (Al+O),  (A4l),  and  (A42)  the  following  expres¬ 
sions  for  and  may  be  obtained: 


(p)  ^p-1  /inrt\r  .  h  /mrt\r  n  Sech  ' 

%  - -V  (f)?p  tMh  (t)Pp  +  +  tMh  («)Pi 

(f)pP 


(p  =  1,2,  .  .  .  r)(b0  »  a)  (M3) 


(r+l) 


br+l 


(l  -  0N)  (^)Pr+l  tanh  (20Pr+l  +  ^JT+1^ 


(r  »  l)  (M4) 


The  rate  of  convergence  of  the  serieB  of  frequency  equations  given 
in  equation  (A3 7)  is  increased  by  substr acting  the  expression 


4  2 
2%  *6 


and  adding  the  equivalent  closed-form  expression 

■  ■) 

The  resulting  equation  is 

kB2  tan  kBkg  +  2kB3kg5  ^  - - - -i - 

°-X’3'5  **  §  %(l)  - 


M* J 


0  (M 5) 
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where  ^  is  determined  by  equations  (M3)  and  (A44) .  It  should  be 
noted  that  equation  (All-5)  will  hold  for  N  stringers  where  N  /  0. 

For  the  case  when  N  ■  0,  handling  equations  (A12),  (A13),  and  (A10) 
when  *  0  in  a  manner  similar  to  that  used  for  the  case  of 

N  stringers  yields  equation  (M5)  except  that 


(M6) 


The  values  P]_  and  are  given  in  equations  (Al4b)  and  (AlJic), 

respectively;  the  value  of  b^  is,  of  course,  b. 


In  order  that  equation  (M5)  together  with  (M6)  may  apply  to  the 
case  for  N  *  ®,  P^  may  be  redefined  by 


B>V 

GL2tc 


(mt) 


Antisymmetrical  Transverse  Modes  of  a  Free-Free  Multistringer  Beam 

Appropriate  trigonometric  series  for  the  antisymmetrical  transverse 
modes  of  a  free-free  beam  are 


up(x,yp)  »  amn(p)  C0B  §7  cos  ^ 

m-0,2,4  n-0,1,2  2L  °P 

(p  *  1,2,  .  .  .  r  +  l)  (A48) 


nuoc 


w(x)  -  Cx  +  >  cm  sin  ^ 

m*2 ,4 ,6  J 


(M9) 


As  in  the  case  for  the  symmetrical  modes  the  choice  of  the  particular 
series  was  guided  by  the  orthogonality  required  for  simplification  of 
the  expressions  in  the  strain-energy  equation;  the  term  Gx  was  added 
in  order  to  allow  sufficient  freedom  of  w(iL). 


.  -  .•«&’ 
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Treatment  similar  to  that  accorded  the  case  of  the  symmetrical  modes 
yields  the  equation 


kB 


1  -  kgkg  cot  kgkg  +  2kg  kg 


4,.  4  £ 


”■=2.4,6  ^2 


(1) . 


*ff)‘ 


0  (A50) 


where  kg,  kg,  and  k^  are  given  in  equations  (Al4f),  (Al4g),  and  (Alifj), 

respectively,  and  is  given  hy  equations  (A43)  and  (A44).  As 

before,  this  equation  holds  for  N  stringers,  where  N  /  0.  When  N  *  0, 

equation  (A50)  may  be  used,  where  is  defined  by  (A46)j  for  N  *  », 

the  same  equations  are  used  with  the  definition  of  as  given  in 
equation  (A47 ) . 
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TABLE  I 

EFFECT  OF  STRINGER  DISTRIBUTION  ON  FREQUENCY' 


N  =  4 


N  =  oo 


Solution 

*2 

1st  symm. 
mode 

2d  symm. 
mode 

3d  symm. 
mode 

Elementary 

+ 

transverse  shear 

5-42 

24.5 

48.3 

Substitute- 
stringer,  bg/bg  *  0*5 

5.16 

20.1 

37-8 

Multistringer 

5- 07 

19.7 

38.1 

Elementary 

+ 

transverse  shear 

5.42 

24.5 

48.3 

Substitute- 
stringer,  bg/bg  =0.5 

5-22 

20.8 

39-1 

Multistringer 

5-15 

20.2 

38.8 

Elementary 

+ 

transverse  shear 

5.42 

24-3 

48.3 

Substitute- 
stringer,  bg/bg  =0.5 

5-22 

20.9 

39.2 

Multistringer 

5.16 

20.4 

39-1 

Elementary 

+ 

transverse  shear 

5-42 

24.5 

48.3 

Substitute- 
stringer,  bg/bg  *  0.5 

5.23 

21.0 

39.5 

Multistringer 

5.19 

20.8 

39.9 

A|«  A2  ■  A3  •  A4* 

b|»b2»b3«b4  «|- 


tc  bfc 


(a)  Multistringer  beams. 


n 


Ap«Ao,  AL«A+btC 

(b)  Substitute-stringer  beam. 


Figure  2.-  Cross  sections  of  multistringer  beams  and  substitute-stringer 

beam. 


Figure  J.-  Effect  of  substitute-stringer  location  on  the  accuracy  of  the 
substitute -stringer  approach  for  box  beams  with  different  stringer 
distributions.  The  labels  IS,  23,  and  5S  denote  symmetrical  modes. 


Cot*  II . 


Figure 


4.-  Cross  section  of  doubly  symmetrical  multistringer  box  beam 
for  either  even  or  odd  number  of  stringers  in  a  cover. 
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